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Abstract 

We compute the one-loop anomalous dimension for the light cone 
distribution function of a heavy quark and solve the corresponding 
evolution equation analytically. Some implications of the results for 
inclusive B decays are discussed. 



1 Introduction 



The heavy mass expansion for inclusive decays of heavy hadrons has become 
a generally accepted tool. Although there are hints that there could be prob- 
lems in the calculation of lifetimes (the A;,-lifetime) within this framework, 
the l/mq expansion seems to be an appropriate method to access semilep- 
tonic processes [|l], 0, 0. 

However, if one is aiming at a complete description of the spectra (i.e. 
energy of the charged lepton or the hadronic invariant mass \/Q^) in 
general problems are encountered in the endpoint region, where becomes 
small. Typically this invariant mass is of the order rriB, the mass of the de- 
caying hadron. This is true in almost all phase space; however, getting closer 
to the endpoint it has become popular to distinguish two regions. The first 
region is the one with very small Q^^j ~ ^qcd also very small hadronic 
energy v ■ Qh ~ ^qcd (here v is the velocity of the decaying hadron); here 
the proper description is a summation over the very few exclusive channels 
contributing in this region. 

However, in the second region where the hadronic energy is still large 
v ■ Qh ~ but the invariant mass becomes small, namely ~ 
^QCD^Bi there is still the possibility to obtain a proper QCD description of 
the decay. In this case it has been shown |^, ^ that one may resum the 
leading twist terms into a universal light-cone distribution function (or shape 
function) for the heavy quark. This function is entirely non-perturbative and 
may be written formally as a forward matrix element of a non-local operator. 
This nonlocal operator corresponds to the Fourier transform of the Wilson 
line ii. 

On the other hand, there are computable perturbative corrections which 
exhibit logarithmic singularities close to the endpoint region. To first order 
in Us the leading contribution to the spectrum close to the endpoint is of the 
form asln(Q^)/Q^, where is now the partonic analogue of the hadronic 
invariant mass of the final state. Upon integration this yields the usual 
Sudakov logarithms. 

Several methods have been proposed to combine both perturbative and 
non-perturbative contributions. At least for the doubly logarithmic terms 
this should be possible since they are process independent and hence could 
be assigned to the universal shape function. In addition, close to the end- 
point, a resummation of the double logarithms becomes necessary, since these 
terms become large. Such a resummation may be achieved by setting up an 
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evolution equation for the shape function as it has been proposed in [|, |^ . 

In the present note we use such an evolution equation and solve it an- 
alytically by putting in the one loop result for the anomalous dimensions. 
In the next section we express the hadronic invariant mass spectra in terms 
of the shape function, thereby defining a limit in which the leading twist 
dominates. In section 3 we consider the renormalization of the shape func- 
tion and set up an evolution equation and derive the anomalous dimension 
in section 4. Finally we solve the evolution equation analytically and discuss 
the implications for the evolution of the moments of the shape function. 

2 Large energy limit of hadronic invariant 
mass spectra 

In inclusive decays such as -B ^ + 7 and B — > Xu + I + ui the effective 
hamiltonian takes the general form 

neff = J,H^, (1) 

where if^ is a hadronic current and is either leptonic or photonic. 

We shall in the following consider the invariant mass spectrum of the final 
state hadrons. We denote with Q\ the square of this invariant mass. The 
spectrum may be written as 

^ = ^ / c^<f . E / d^x{27:Y6\P - q - Px)KQl - PI)L,^H^^, 

(2) 

where phase space integration over the leptonic and hadronic total momen- 
tum Px is denoted by respectively. 

The leptonic (or photonic) and hadronic tensors are defined as 

V = (3) 
H,, = {B{P)\H,\X){X\Hl\B{P)). (4) 

The non-hadronic part can be calculated perturbatively and is treated in the 
following as a known function of q. The object of our interest is the hadronic 
part, 

Qf.u = J2 f d^x{27r)^S\P - q - Px)H^,, (5) 
X 
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which in general cannot be determined perturbatively. Proceeding along 
standard lines we have 

g^. = J d'xe'^^{B{P)\Hl{0)H,{x)\B{P)) (6) 

The hadronic current is usually a bilinear function of two quark operators 
one of which is a heavy b quark. Using the fact that the heavy mass sets a 
large scale compared to Aqcd we may set up an operator product expansion 
in the usual way. We write 

= q{x)r^b{x) = e-'"'^'^q{x)rM^), (7) 

where q{x) represents a massless quark and the large part of the 6-quark mo- 
mentum has been explicitely removed from the field by a phase redefinition. 
This leaves us with 

g^.- J ci^xe-^^("^^''-^)(5(P)|6,(0)rtg(0)g(x)rA(x)|i^(P)), (8) 

In the following we are interested in the endpoint region, which is defined 
by a specific kinematical limit denoted large-energy limit. We shall define 
this limit first in terms of the partonic total momentum of the final state 

= mbyf" - (9) 

corresponding to the partonic invariant mass and the partonic energy 
V ■ Q. In the endpoint region the total energy v ■ Q of the final state scales 
with the heavy quark mass nib while the invariant mass is of 0{AQCD'nT'b) 
such that the light-cone component k+ of the total final state momentum 
remains finite of 0{Aqcd)' 

nib 

v-Q 

2vQ 

nib 



oo 



const. (10) 



We shall relate the partonic variables to the hadronic ones at the end of the 
section. 

It has been shown that in the limit ( [T0| ) the leading contribution may 
be obtained by contracting the light quark 

(5(P)|6,(o)rt^r,Pexp 

This expression may be Taylor-expanded around x = and the expansion 
resummed to yield formally 

Q.. = / ^e(fco)5(fc^) 

(S(P)|6,(0)rt ^r,(27r)454(g -k + iD)hM\B{P)) ■ (12) 

In the following, 5-functions as in the preceding expression are defined by 
their Fourier transform, which involves a path-ordered exponential as in (|1T]) . 
Performing the fc-integration we get 

Q^, = (5(P)|6,(0)rt((^ + z^)r,5((Q + zD)2)6,(0)|5(P)). (13) 

In the endpoint region (p!OD Q is almost light-like and can be decomposed 

as 

Q = {v- Q)n+ - k with vn+ = l, nl = 0, k = C(Aqcd) • (14) 
Thus in (|13|) we can approximate (i-D+ = n+ ■ iD) 

{Q^iDf ^ Q^ + 2vQiD+ 

h rti h 

neglecting terms of O^KqcD/^b)'- 

Q^, = {B{v)\Kml(^^J{Q^ + ■ Q^D+)KmB{v)) (15) 

For the perturbative calculation in the next section it is useful to rewrite Q^,^ 
in terms of the imaginary part 

Q^, = Im-f^, , (16) 



X 

-I j dx ■ A{x) 

'0 



6.(a;)|P(P)). (11) 
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where 

f„„ = -(B(o)\K(«K<S!T, ^^ ^ '^^^^ ^ ^ KmB(v)) . (17) 

Using heavy-quark symmetry to disentangle the Dirac structure of this 
expression we finally arrive at 

g^. = lTr{rt^r,p+} j dk+diQ' + 2v ■ Qk+)f{k+), (is) 

where Py = (1 + ^)/2. 

Here we have introduced the shape function or light-cone distribution 
function for the heavy quark 

/(fc+) = {B{v)\hJ{k+ - iD+)K\B{v)) . (19) 

The light-cone structure function f{k^) measures the probability of finding 
a heavy quark with light-cone component k^ inside the B~meson and is a 
universal function. 

Up to now we have worked in terms of partonic variables, while exper- 
imentally only hadronic variables are of interest. Introducing the hadronic 
light-cone variable in terms of the hadronic invariant mass squared and 
the hadronic energy v ■ Qh we have 

+ 2vQh 2{vQ)+2A + K / bj K J 

wc find that the support of f{k+) is the interval — oo < k^ < A such that 
< 0. In reality, the spectrum is peaked near k^ = 0. The shape function 
may directly be measured in a inclusive semileptonic or radiative decay by 
measuring the spectrum of the hadronic light-cone variable K^. 



3 Renormalization of the light— cone structure 
function 

In the last section we have shown, that in the large-energy hmit the inte- 
grated hadronic tensor Q^i, can be approximated by the quantity Q/j^i, defined 
in the large energy limit of QCD. Up to now only non-perturbative correc- 
tions resummed in the structure function are included. 
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To really establish the approximation we have to check that this also 
works perturbatively, since it is well known that in the kinematical region 
of small invariant mass the perturbative spectrum is plagued by large 
logarithmic corrections. Therefore we should require that order by order in 
perturbation theory the leading IR singularities of Q^jj arising in the large- 
energy limit should be reproduced by Q^u- The leading logarithms are uni- 
versal and hence we may assign them to be perturbative corrections to the 
structure function. 

In order to identify the relevant IR singularities which have to be re- 
produced by perturbative corrections to the structure function we shall first 
consider the hadronic invariant mass spectrum calculated in full QCD in the 
large-energy limit. The relevant terms of the one loop contribution can be 
generically written as 



^r(per-i)(l) ^ 



+ 



(21) 



where the real functions fi{p) are regular in the limit p ^ with A;+ fixed. 
From this we easily identify the IR singularity as the ln(Q2/m^)-term since 
k+ is held fixed in the large-energy limit. If the large-energy limit makes 
sense in perturbation theory, we have to require that exactly this term is 
reproduced by the (9(a<i)-corrections to the light-cone distribution function. 

The origin of this logarithmic divergence in the spectrum can be traced 
back to double logarithms of the correlator T^^^ . Thus to one- loop order the 
double logarithmic terms in T^^ and T^^ have to match in the large-energy 
limit. It therefore suffices to look for the most singular Feynman integrals 
contributing to the correlators at 0{as)- 

The one-loop corrections to the shape function consist of vertex diagrams, 
self-energy contributions and a box-type diagram. However, in full QCD the 
UV divergencies of the self-energy diagrams cancel with the UV divergency of 
the vertex diagram because of current conservation. For the shape function 
divergent pieces appear in the self energy of the heavy quark which will 
be taken into account below. The self-energy insertion into the light-cone 
propagator does not contribute. Thus, we can restrict ourselves to the vertex 
diagrams shown in figure ^ 
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Figure 1: Renormalization of the structure function. The plain line denotes 
a massive or massless QCD propagator, the double plain line a heavy quark 
propagator, the dotted line a light-cone propagator. The shaded blob on 
the r.h.s. symbolizes the insertion of the structure function as a nonlocal 
operator. 



The part of the vertex correction in full QCD (first term in Fig.|l|) respon- 
sible for double logarithms reads 

- (22) 



d 



Di 



(27r)« iv-l + vk' + + k'Y)l^ [(£ + k' + QY + le] ' 

where k' is the residual momentum of the heavy quark and the integration 
measure indicates dimensional regularization with D = 4 — e. 

On the other hand, we can directly compute the corrections to the shape 
function either by reading off the propagator from (0) or by the usual way 
of calculating corrections to Wilson lines 0. The corresponding expression 
is 



d 



(27r)^ {vt + v k')P{Q^ + 2v ■ Q{i+ + k'+) + ie) 



in which the leading logarithms have to match the full QCD result (|2^). 
Evaluating the corresponding integrals / and / to double logarithmic 
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accuracy and supressing constant terms we get 

' 1 _\,^^.fK^K^\ (24) 



' (4^)^-^r(i-^)r(e)r(i)f '-^-'^-" ^) (25) 



(47r)2' ' v-Q ^ ' ' ' 

1 [2 2 f — k', — ie\ c r, f k+ — k'^ — ie 
- In H h m 



(47r)2 V ■ Q 



e^ey yu ye \ 



where 



ln(47r) - 7s . (26) 



Note that is vahd up to single logarithmic corrections. 

Due to the presence of both UV and collinear divergencies, a pole 
appears. However, both singularities are of UV nature in the sense that the 
integral converges for e > and can therefore be removed by renormalization, 
as will be explained below. 

In the difference of both integrals 

AI = I-i (27) 

2 2, f k< — k\ — ie\ c 1 , r, f k+ — k'_j_ — ie 
' In — In^ 



(47r)2 V ■ Q 



s"^ ' e \ fi J e 2 \ 

the double logarithms do not cancel. 

The operator to be renormalized [i.e., the operator defining the light-cone 
distribution function ([T9|)] is of nonlocal nature. Hence, we can introduce an 
integration over k'_^_, which plays the role of a summation over a continuous 
operator basis labeled by k'_^_, to absorb the local as well as the nonlocal UV 



divergencies in (p7|) into a renormalized distribution function. 

Furthermore, the difference of integrals ( JI7\ ) has a finite part. In order 
to reproduce the amplitudes (defined in a definite scheme, e.g. MS) of full 
QCD in the effective theory, these parts, the matching corrections, have to 
be included in the renormalization constants. In ordinary effective theories 
such as HQET, such finite matching corrections will not affect the scaling 
behavior to one-loop order, so they have to be taken into account explicitly 
only from two- loop order on. By contrast, in the present case the finite 
matching corrections 

k4- — k': — ie^ 
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X / dk'. 



explicitly depend on the renormalization scale jj, and therefore modify the 
scaling behavior already at one-loop order[|. 

In order to find the renormalization kernel, we take the imaginary part of 
the difference of vertex corrections to the correlators T and T resulting from 
([27|), keeping only UV-poles and double logarithms: 

AQ^^ = ^Im(l^(rng) - (29) 
= ^^Tr{r^t^rP+} /rfM(Q' + 2t;-gA;+) (30) 

TT 2 J 

-Z^'\k+,k'^,fi)MS' - —KK - k'Af^'^'-'Kk'^) 

t J 

The second term in square brackets accounts for wave function renormaliza- 
tion of the heavy-quark fields. The one-loop contribution to vertex renor- 
malization is given by Q 

-c.^A[e(.;-Min(^)] 

+ ^C^i^^[e(fc;-Mln^(^^^)], (31) 

which has to be understood in the distribution sense. 

In this way we can define a renormalized shape function 

/(fc+,/i) = y d}^^Z{k^,]^^,^)MS'f'"'''\K) (32) 
where the one loop renormalization kernel reads 

Z{k+,k'^,^i)MS' = - + {-^)Z^^\k^,k'^,^i)MS'- (33) 



^The appeareance of nonlocal divergencies and double logarithms in the matching cor- 
rections spoils the naive formulation of a large-energy effective theory (LEET) for exclusive 
processes, as has been observed in |l2| . 

^The renormalization prescription which matches QCD in the MS scheme is denoted 
as MS'. 
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In what follows it will be convenient to switch to the MS''-scheme setting 

l = l-Un{A7i) + ]pE (34) 
and keeping terms proportional to and 1/e in the renormalizat ion-kernel: 

- c^F-;7^[e(fc;-MM^i^)] 

e dkj^ /i 
+ ^^^Jr[0(fc'.-Mln^(^^^)] (35) 

Then to one loop order 

Q^, = iTr{ri^r,P+} j dk+6{Q^ + 2v ■ Qk+)f{k+, M) (36) 

is UV-finite. Here M is a scale of the order v ■ Q which is a perturbative 
scale. As a consequence of the finite renormalization in ( pSj) also the correct 
collinear singularity is reproduced. 

4 Evolution equation 

As a consequence of the renormalization procedure the structure function 
now depends on both the renormalization scale /x and the renormalization- 
scheme {MS ) and obeys the evolution equation 

^/(A;+,/i) = j rfA;;r(A;+,fc;,/i)/(fc;,/i). (37) 

The evolution kernel T{k+,k'_^_, /j) in the MS''-scheme0 is defined implicitly 
by 

J dklT{k+, k'l, fi)Z{kl, k'^, /x)^' = ^^]^^(^+' 1^)ms'^ (38) 
■^In what follows this scheme dependence of the evolution kernel is understood. 
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where 



« .2a.r^ff„ + ...)^ (39) 



dlnfi (91n/i Vtt / 9a 

with Pq = (33 — 2nj)/12. Thus, to one-loop order the evolution kernel is 
given by the negative coefficient of the 1/e-pole term in Z^^\k+, k'_^_, jj,)-^' 
and the partial ln(/i)-derivative of its finite partQ 

r(A;+,A;;,/i) = (^)r«(A;+, fc^, /i) (41) 

TT 

= i^)c, (^sik'^ - k^) + ^m'+ - K) 

In order to solve the evolution equation we Fourier transform both the struc- 
ture function and the evolution kernel: 

f«(e,/.) = |^r«(A:+,0,/.)e'^-^ (42) 

= ^{i-in(ie/i|)-^f(e(0-e(-0)} (43) 

Then the evolution equation (^) reads 

/(e,/^) = (-)(2vr)f«(e,/x)/(e,/^) (44) 



din fx TT 
and can be solved easily 

h^,f^) = Ui^,l^,l^o)m,^o) (45) 

where 

U (e, /i, /io) = K{i2, fio) exp (i^^^iO) (46) 



'*Note that the double pole in k'_^, lJ')jjg' cancels with the explicit /i-derivative 



of the single pole-term: 



_l__^[e(,;-, J ,.„^i_i±l,|. (40) 
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with 



u; = ^iJ^l^], e(0 = e(e)-e(-o (47) 



2/3o \as{no 
and 

=(1LY^(^Y^.^'^V^^\ (48) 

V/io/ V«4/^o)/ 

Transforming back to momentum space we get 



/i) = y dk'_,u{k+, k'^- /i, /io)/(fc;, /io) (49) 

where 

f/(fc+, fc;; /i, /.o) = -i^'(/i, /xo)r(i + ^k ''^'' (fcf-fcji+^ • 

Note that the tree-level result is recovered from (^U|) by letting /i = /iq and 
performing the limit uj —0 carefully, since 

Mm -'---"^(K-K) ,^), (51) 

Once radiative corrections in the definition of the light-cone structure 
function are included, we are in the position to study the invariant mass 
spectrum. To this end we insert (^) into (PB|), choose for /i a high energy 
scale M = 0{y ■ Q) and let = fj, 

Q,, = iTr{rt ^r,P+} J dk'^Uik^, k'^- M, /i)/(A;;, /i) (52) 

The coefficient function U{k^, k'^; M, fi) is analogous to the set of Wilson 
coefficients Ci{M,fi) in local effective field theories. 

Note that, by definition, the convolution of the renormalized structure 
function with the coefficient function is independent of the choice of the 
renormalization scale fi. At the low energy scale = 0{Aqcd) the structure 
function is only weakly affected by perturbative corrections and can be safely 
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considered as a pure non-perturbative object which may be described by a 
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f{k+,AQCD) = f^'^'"'''\k^) (53) 

Putting everything together, the renormahzation group improved hadronic 
tensor becomes: 

= -^pi-^/f(Af,,»)r(i + i^)f."™TV{rt«r„p+} (54) 

(88) 

5 Evolution of Moments 

It is well known that the moments of the light-cone structure function may 
be related to matrix elements of local operators. The moments satisfy the 
relation 

Mn = j dk+{k+rf{k+) = {B{v)\K{iD+rK\B{v)) (56) 

in particular we have Mq = 1 and Mi = due to the normalization and the 
equations of motion. 

However, since the light-cone distribution funtion evolves with a change 
of scale, also the moments depend on the scale. The moments at the scale 
are given by 



MM = J dk+{k+rf{k+,fi) = J dk+ik+r J dk'_^u{k^,k'^,fi,fio)f{k'+,t^o) 

(57) 

Since the kernel U{k+, k'_^_, /i, /io) depends only on the difference of /c+ and k'_^_ 
one may reexpress these moments in terms of the moments at scale /Xq by a 
simple change of variables. One obtains 

Mnifi) = V ..M fio) / dk+{k+r-^U{k+,0,fi,fio) (58) 

which means that under renormahzation the n*^ moment will depend on all 
the moments Mj with j < n at the lower scale. 
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However, due to the power type behaviour in the variable — k'^ of the 
evolution kernel U{k+, k'_^, /i, /iq) the integral in (^) does not converge, since 
the region of integration is — cxd < A;+ < and hence it has to be regularized. 
If we regularize the integral by restricting the integration region to —A < 
A;+ < with some cut off A, the dependence on A can be already guessed 
from dimensional analysis. More precisely, the moments M„(/i) exhibit power 
divergencies where the dependence on A is given by 

n 
j=0 

where from (pOD we get 

^ , , n\ sinvrcij 

J- ['^ J )■ ^ 

To interpret this result it is instructive to consider the first few moments. 
The normalization receives a multiplicative renormalization 

Mo(/i) = Mo(/io)Co,o(/i,/io) [j-Y (61) 
while the first moment receives contributions from both Mo(/io) and Mi(/io): 
Mi(^) = Mi(/io)Ci,i(/i,/io) (^)" + AMo(/io)Ci,o(/i,/io) (y)" (62) 
However, the first moment is 

Ml = {B{v)\K{iv ■ D)K\B{v)) (63) 

which according to the equations of motion should vanish. This is true, if 
the pole mass of the heavy quark is taken as the expansion parameter, such 
that no residual mass term appears. The choice of a different mass definition 
m ^ m + 6m will modify the equation of motion into 

{iv ■ D)hy = 5mhy (64) 

and hence the running of the first moment implies a change in the pole mass 
definition due to the power divergence 

<5m = ACi(/i,/io) . (65) 
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In principle this relation can be used to fix the value of the cut-off in terms 
of the pole mass; once this is done, all higher moments can be computed in 
terms of the pole mass at the lower scale fiQ. 

Although the shape function is an object entirely defined in HQET and 
hence should be independent of the mass m^, in this indirect way a mass 
dependence comes into the game, at least if we insist to interpret the moments 
in terms of the matrix elements (|56D . 

6 Conclusions 

In this paper we have proposed a method to combine perturbative and non- 
perturbative contributions to hadronic invariant mass spectra in the endpoint 
region where the invariant mass becomes small of the order 0{Aqcd^b)- 
At tree level the non-perturbative corrections of leading twist are resummed 
into a universal light-cone distribution function. However, taking radiative 
corrections into account, this function has to be renormalized. In order to 
match the leading IR-singularity of the QCD spectrum, we have chosen a 
renormalization scheme which except standard UV-renormalization provides 
corrections of the IR-behaviour of the structure function adding proper finite 
terms to the renormalization kernel. As a consequence of renormalization the 
renormalization scale dependence of the structure function is described by 
an evolution equation. We have computed the evolution kernel to one loop 
order, which deviates from the usual MS'-kernel, since we have included a 
finite renormalization in order to reproduce the infrared behaviour of full 
QCD. 

The analytical solution of the evolution equation yields a resummation 
of logarithmic corrections and relates structure functions at different renor- 
malization scales. The last property has then be used to include radiative 
corrections in the invariant mass spectrum in a manner common from usual 
renormalization group techniques. 

As the structure function itself is scale dependent so are its moments. 
However, including radiative corrections the moments exhibit an UV-diver- 
gence which has been regularized by an hard cut off thereby mixing moments 
of different order. We fixed the cut off by relating the first moment to the 
definition of the pole mass. Once the value of the cut off is known, the scale 
dependence of all other moments is computable. 



15 



Acknowledgements 

We are grateful to M. Neubert and V. Smirnov for valuable discussions. This 
work was supported by the "Graduiertenkolleg: Elementarteilchenphysik and 
Beschleunigern" and the "Forschergruppe: Quantenfeldtheorie, Computeral- 
gebra und Monte Carlo Simulationen" of the Deutsche Forschungsgemein- 
schaft. W.K. is supported by German Bundesministerium fiir Bildung und 
Forschung (BMBF), Contract No. 05 6HD 91 P(0). 

References 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 

[7 
[8 
[9 

[lo: 
[11 

[12 
[13 



A.F. Falk, M. Luke and M.J. Savage, Phys. Rev. D53, 2491 (1996). 

A.F. Falk, M. Luke and M.J. Savage, Phys. Rev. D53, 6316 (1996). 

A.F. Falk and M. Luke, Phys. Rev. D57, 424 (1998) 

M. Neubert, Phys. Rev. D49, 3392 and 4623 (1994). 

T. Manuel and M. Neubert, Phys. Rev. D50, 2037 (1994). 

I.I. Bigi, M.A. Shifman, N.G. Uraltsev and A.I. Vainshtein, Phys. Lett. 
B328, 431 (1994) 

G.P. Korchemsky and G. Sterman, Phys. Lett. B340, 96 (1994). 
A.G. Grozin and M. Neubert, Phys. Rev. D55, 272 (1997). 
A.G. Grozin and G.P. Korchemsky, Phys. Rev. D53, 1378 (1996). 
M.J. Dugan and B. Grinstein, Phys.Lett. B255, 583 (1991). 
U. Aglietti and G. Corbo, |hep-ph 9803485| . 



C. Balzereit, T. Manuel and V. Smirnov, in preparation. 
The BABAR Physics Book, SLAC-R-504, in preparation. 



16 



